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Hennebelle & Chabrier ( 2008 ) (HC08) attempted to derive the stellar IMF as a consequence of log- 
normal density fluctuations in a turbulent medium, using an argument similar to |Press & Schechter ( 1974) 
for Gaussian random fields. Like that example, however, the solution there does not resolve the "cloud in 
cloud" problem; it also does not extend to the large scales that dominate the velocity and density fluctua- 
tions. In principle, these can change the results at the order-of-magnitude level or more. In this paper, we 
use the results from Hopkins (201 1 ) (HI 1) to generalize the excursion set formalism and derive the exact 
solution in this regime. We argue that the stellar IMF and core mass function (CMF) should be associated 
with the last-crossing distribution, i.e. the mass spectrum of bound objects defined on the smallest scale 
on which they are self-gravitating. This differs from the first-crossing distribution (mass function on the 
largest self-gravitating scale) which is defined in cosmological applications and which Hll show corre- 
sponds to the GMC mass function in disks. We derive an analytic equation for the last-crossing distribution 
that can be applied for an arbitrary collapse threshold shape in ISM and cosmological studies. With this, 
we show that the same model that predicts the GMC mass function and large-scale structure of galaxy 
disks also predicts the CMF - and by extrapolation stellar IMF - in good agreement with observations. 
The only adjustable parameter in the model is the turbulent velocity power spectrum, which in the range 
p w 5/3 — 2 gives similar results. We also use this to formally justify why the approximate solution in 
HC08 is reasonable (up to a normalization constant) over the mass range of the CMF/IMF; however there 
are significant corrections at intermediate and high masses. We discuss how the exact solutions here can 
be used to predict additional quantities such as the clustering of stars, and embedded into time-dependent 
models that follow density fluctuations, fragmentation, mergers, and successive generations of star forma- 
tion. 
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1 INTRODUCTION 

The origin of the stellar initial mass function (IMF) is a question 
of fundamental importance for the study of star formation, stellar 
evolution and feedback, and galaxy formation. It is an input into a 
huge range of models of all of these phenomena, and a necessary 
assumption when deriving physical parameters from many obser- 
vations. However, despite decades of theoretical study, it remains 
poorly understood. A critical first step - although by no means a 
complete description of the origin of the IMF - is understanding 
the origin and form of the mass function of protostellar cores (the 
CMF), specifically that of self-gravitating, collapsing cores that 
will ultimately form stars. 

Recently, HC08 presented a compelling argument for the 
physical origin of the IMF shape, as a consequence of the CMF re- 
sulting from lognormal density fluctuations in a turbulent medium. 
It is increasingly clear that the density structure of the ISM is dom- 
inated by supersonic turbulence over a wide range of scales (e.g. 
Scalo & Elmegr een|2004[|Elmegreen & Scalo|2004l|Mac Low &| 
Kles sen|2004||McKee & Ostriker|2007| >, and a fairly generic con- 



this and cosmological Gaussian density fluctuations, HC08, build- 
ing on the earlier work in|Inutsuka l |2001) as well as |Padoan et al.| 
l |1997) ; |Padoan & Nordlund| ( |2002 i, attempted to approximate the 
mass function of self-gravitating cores in a manner exactly analo- 
gous to Press & Schechter| ( |1974) . If the density field <5(x) (where 
8 oc lnp) at the point x is normally distributed, then the average 8 
smoothed around the point x with the appropriate window function 
of effective radius R, 8(x\R) is also normally distributed, with a 
variance S — cr~(R) that can (in principle) be calculated from the 
power spectrum or simply estimated from the Mach number (see 



e.g. Passot & Vazquez-Semadeni[1998| 


Nordlund & Padoan 1999; 


Federrath et al.||2008| |Price et al.|2011 


i. The total mass which is 



sequence of this is that the density distribution converges towards 
a lognormal PDF, with a dispersion that scales weakly with Mach 
number (e. g. |Vaz quez-Semadeni 1994; Pad oan et al.|1997| |Scalo| 
|et al.|[T998) |Ostriker et al.||199*9*| . Based on the analogy between 
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self-gravitating, over the scale R, is then just oc f^L p(S) P(S)dS, 
where B(R) is the critical density above which the region would be 
self-gravitating. Differentiating this total mass fraction with respect 
to the mass scale associated with each R gives - approximately - 
the total mass in bound objects per unit mass, M 2 dN /dM, hence 
the mass function. HC08 showed that, for plausible 5, and B given 
by the Jeans condition for thermal plus turbulent velocities, this 
argument reproduces all of the key features of the observed CMF 
and stellar IMF (for extensions of this calculation allowing for dif- 
ferent gas properties and , see Hennebelle & Chabrier 2009, |20lTj 
IChabrier & Hennebelle|20TTT l. 

While extremely interesting, there are, however, a number of 
uncertain assumptions in this derivation of the CMF. Because it 
focuses on small scales exclusively, the properties of the "parent" 
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scales (e.g. GMCs) must be assumed somewhat ad hoc. Since most 
of the power in velocity (hence density) fluctuations is on large 
scales (true for any reasonable turbulent power spectrum, and also 
observed; see Ossenkopf & Mac Low 2002| |Brunt et al.|2009| >, this 
means S(R) could also not be calculated but was instead assumed, 
and its "run" with radius R (being undetermined) was neglected 
(and the resulting mass function artificially truncated before going 
to very large scales). 

Most important, like Press & Schechter ( 1974), this derivation 
does not resolve the "cloud in cloud" problem. A given region may 
well be self-gravitating on many different scales R (the smoothed 
5(x\R) crossing back and forth across the critical B(R) with scale). 
This makes the resulting mass function inherently ambiguous, since 
different crossings of the same spatial location are counted multi- 
ple times and with different signs. [Bond et al.| ( |1991| > resolved this 
ambiguity by extending the mathematical excursion set formalism 
and defining the "first crossing distribution." This allowed them to 
rigorously calculate the mass function of objects - counting each 
point only once - where "mass" was defined as the mass enclosed 
in the largest scale R on which a region was self-gravitating. 

More recently, HI 1 showed how the full excursion set formal- 
ism could be generalized to the problem of lognormal density fluc- 
tuations in a gaseous galactic disk. By including disk scale-effects, 
they showed that it is also possible to predict the absolute mass 
scales, variance a~(R), and barrier B(R) - quantities which had to 
be assumed ad hoc in HC08 - with only the assumption of a turbu- 
lent spectral shape. On small scales (below the disk scale height), 
the collapse condition is just the Jeans condition for turbulent plus 
thermal support - identical to HC08- so there is no reason why the 
approach therein cannot be extended to the same scales. 

Hll showed that the first-crossing distribution for turbulent 
gas in a galaxy disk is not the CMF or stellar IMF, but rather agrees 
extremely well with observations of the mass function and other 
properties of giant molecular clouds (GMCs). This highlights just 
how important the distinction of multiple crossings can be - ex- 
tending the HC08 argument to large scales, self-gravitating objects 
with scales ~ 10 6 Mq would be recovered (and contain more mass 
than the objects self-gravitating on scales ~ IMq). Clearly, these 
are the "parent" clouds, not the protostellar cores! 

So how can the two be distinguished? Physically, consider 
a region which is self-gravitating on a large scale Rq. If it con- 
tains multiple sub-regions that are themselves self-gravitating on 
a smaller scale Ri, then the entire Ro object will not form a sin- 
gle core (see e.g. the discussion of the "cloud in cloud" problem 
in |Veltchev et al. 201 1 ; Do nkov et al.| 201 1 1. Since the mean den- 
sity at Ri to be self-gravitating must be larger than that at Ro, these 
sub-regions will collapse more rapidly - the "parent cloud" is frag- 
menting into smaller objects. This can be continued iteratively in- 
side Ri. It is only when a region is self-gravitating on a scale Ri, 
and not self-gravitating on any smaller scales, that its collapse will 
proceed without fragmentation. 

Therefore, we argue in this paper that the CMF (and, to the 
extent that it is related, the IMF) should be associated with the last- 
downcrossing distribution: i.e. the mass function of regions which 
are self-gravitating, but with mass defined at the smallest scales 
on which they are self-gravitating. This has not, to our knowledge, 
been studied in any cosmological context (since halos are assumed 
to collapse "into" the mass they contain), so in § [2] we derive a 
rigorous analytic expression for this mass function as a function 
of arbitrary collapse threshold. In § [3] we combine this with the 
model from Hll, which gives the appropriate collapse threshold 
and variance for a galactic disk. In §[4] we compare the results to 
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Figure 1. Predicted last-crossing mass function (distribution of bound 
masses measured on the smallest scale on which they are self-gravitating) 
for a galactic disk with turbulent spectral index p = 2, and Mach number at 
scale ~ h of Mi, = 30. In units of the sonic mass M son j c = (2/3) c^R son i c /G 
and total disk gas mass M gas , all other properties are completely specified by 
disk stability. We calculate this with the analytic iterative solution to Eq. |14| 
the standard Monte Carlo excursion set method gives an identical result. 
We compare the first-crossing distribution - the distribution of masses mea- 
sured on the largest scale on which gas is self-gravitating, which in HI 1 we 
show agrees extremely well with observed GMC mass functions. The MF 
derived in HC08 by ignoring multiple crossings is also shown. We compare 
the observed stellar IMF from Kroupa 1 2002 1 and Chabrier ( 2003 ) (shifted 
to higher masses by a simple core-to-stellar mass factor of 3). 



the observed CMF and stellar IMF, see how it contrasts with the 
first-crossing distribution (GMC mass function), and examine its 
dependence on turbulent properties. In §[6]we discuss implications 
and future directions for this work. 

2 THE LAST-DOWNCROSSING DISTRIBUTION 

Our derivation here closely follows that of the first-crossing distri- 
bution in Zhang & Hui (2006), to whom we refer for more details. 

Consider the Gaussian field <5(x[i?) (which for this problem 
represents the logarithmic density field smoothed in a kernel of ra- 
dius R about x); the variance in the field = S(R) is a monotonically 
decreasing function of R so we can take S as the independent vari- 
able and consider 5(x \ S). For convenience we will drop the explicit 
notation of x. The PDF of S(S) is just 

p ° (5|5)= ^ exp (-S) (1) 

The barrier B(S) is the minimum value S(S) which defines 
objects of interest (e.g. collapsing regions). Normally, we would 
define the first-crossing distribution by beginning with the initial 
condition <5(0) = at S = (R — > oo), then evaluating S(S) at suc- 
cessively increasing 5 (smaller scales R). This is straightforward: 
given an "initial" value Sq(Sq) at scale So, the probability of a value 
Si(Si) at a scale Si > So is just 

Pw[8i(Si)\5o(So)]=Po(5i-5 \Si-So) . (2) 

This can be integrated until the "trajectory" (or random walk) 5(S) 
first exceeds B(S). 

For the last crossing distribution, we need to determine the 
smallest scale at which a trajectory that had previously crossed the 
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Figure 2. Predicted last-crossing mass function from Fig. [T] compared 
to observed pre-stellar/starless core MFs in different regions. The CMFs 
for Ophiuchus, Perseus, Taurus, and Orion are compiled in Sada voy et al.| 
(2010); we a lso compare the independent determination for Ophiuchus in 
|Enoch et al."]pu08). The CMFs in the Pipe jRathborne et al.|2009) and 
Chameleon I jBelloche et al.|2()l If also agree well but would be indis- 
tinguishable from the other systems on this plot. The normalization of the 
observed/predicted CMF is arbitrary (depending on absolute mass and vol- 
ume) but - unlike the comparison with the stellar IMF - the absolute mass 
value is not rescaled here. Note the dynamic range is limited, and the obser- 
vations have substantial uncertainties, but the agreement is good. 



barrier at some larger scale again falls below B(S). Consider a tra- 
jectory, therefore, in the opposite direction: beginning at an ini- 
tial Sj (Si) at some arbitrarily small Ri — > 0, with a corresponding 
Si(Ri), and evaluating S at successively smaller S (larger R). For 
the last-crossing distribution to be meaningfully defined, it must be 
the case that as R — > 0, the probability of S-, exceeding B(S,) vm- 
ishe&,Pb(B[Si]\Si)->0. 

Define the last crossing distribution fe(S\ <5,)dS as the proba- 
bility that a trajectory, beginning from this state, crosses 6(5) for 
the first time between 5 and S + dS. Also define 11(5 1 St, S) dS as 
the probability for the trajectory beginning at 5i(St) to have a value 
between S and <5 + dS at scale 5, without having crossed B(S) at any 
larger S. Clearly, the integral of ft and integral of II for 5 < B(S) 
must sum to unity: 



1 



B(.V) 



dS' f e (S'\6,)+ / n(<5|<5,-,S)d<5 



(3) 



For a step in this "opposite" direction of decreasing S, the 
probability of going from an initial <5i(Si) to a value 5o(Sa) with 
Si > So is related to the probability of going from So(So) to di(Si) 
by Bayes's theorem: 



Pm [Sq(So) I 5i (Si )] = P l0 [S t (Si ) | <5o(S )] 



Po(So\S ) 



Po(8i\S 1 ) 
P [5 - Si (So/Si) | (Si - So) (So/ St)] 



(4) 



If we ignored the barrier, II(<5|<5,-, S) would be equal to 
Poi [S(S) \ 8i(Si)]. But we must subtract from this the probability that 
a trajectory crosses the barrier at some larger S' > S and then passes 
through S(S), so 

n(5\S„S)=Poi[5(S)\di(Si)]- 



dS' ft (S' I Si) Pot 15(S) 1 5' (S') = B(S')] (5) 



Before going further, note that we do not know the value of 
Sj(Si), but know its distribution. The fraction of trajectories cross- 
ing the barrier at each S is the integral of fe(S\ Si) weighted by the 
distribution of Si. We therefore define 



f e (S) = (f(S\ 5,)) = / dSiP (Si | S,)f e (S\Si) 



(6) 

Technically, the upper limit of this integral should be B(S,), since 
the trajectory must begin as uncollapsed, but since we choose S, 
such that the collapsed fraction is vanishingly small, we can safely 
take B(Si) — > oo. If we take the integral d<5, Po(5; | Si) with re- 
spect to both sides of Equation [5] note that both S and S' are inde- 
pendent of St, and use the fact that 

dSiP (S, | Si) P m [S(S] | Si(S,)] =Po(S\ S) (7) 

we obtain the <5;-averaged equations 

rSi rB(s) 
1 = / dS'f e (S')+ / n(<5|S)d<5 (8) 

J S J — oo 

Il(5\S)=Po(S\S)- [ ' dS'f e (S')P 01 [5(S)\B(S')] (9) 



Taking the derivative of Eq.[8] we obtain 



fi(S) = U(B(S)\S) 



dB 
dS 



B(s) dU(8\S) 
dS 



d<5 



(10) 



Finally, we insert Eq.|9]in Eq.[l0] perform some simplifying 
algebr iQand obtain 

fe(S)= gi (S)+ f S ' dS'f e (S')g 2 (S,S') (14) 
J s 



where 



82 



(s,s')=[ 



(15) 



B(S)-B(S') B(S) „dB 



S-S' 



-2^\P m [B(S)\B(S')] 



(16) 



Equation [14] is a Volterra integral equation, which for a gen- 
eral barrier B(S) has a unique solution that can be calculated by 
standard numerical methods. For example, if we grid S on a mesh 
with equal spacing, 

S„=S,-«AS, n = 0,l,...,N, AS=^^ (17) 
and treat fi(S„) as a vector, then the integral equation becomes a 



1 Specifically, using the relations 

B(S) 



J — OC 



d8P 0l [S(S)\B(S')] 



r B W dP (S\S) 



dS- 



B(S) 



f 

J-oo dS IS 

r B W dp 01 [6(s)\B(s')] 



P (B(S)\S) 



(11) 

(12) 



OS 



A5- 



B(S')-B(S) B(S) 
2(S'-S) 2S 



P 0l [B(S)\B(S')] (13) 
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triangular matrix equation which can be solved iteratively: 

/,(So)=gi(So)=0 (18) 



/*(Si)=gi(Si)(l-fli,i)- 



(19) 



gi(S„) + ft(S n -i)H„,„+Yl [ft(JS m )+ft{S„-i)]H^ 



fe(S„)\„>i = 



where 



AS 
2 



(20) 



(21) 



For a linear barrier, B(S) = Bo + BS, this has a closed-form 
solution: 

fe{S\B = Bo + l3S) = /3P (B(S)\S) = -^L exp f - ^ ) (22) 

Equation [14] is qualitatively similar to the governing equa- 
tion for the first-crossing distribution (compare Eq. 5 in [Zhang 
|& Hui|[2006[ l, but with some critical differences. Up to a sign, 
gi (5) is identical. In the gi term, however, there is an additional 
B(S)/S in the coefficient, and Pqi[B(S) \B(S')] appears instead of 
Pio[B(S) 1 6(5')]. And of course, the integration proceeds in the op- 
posite direction. These corrections result in qualitatively different 
behavior. For example, for the linear barrier, the first-crossing dis- 
tribution has a pre-factor So/5 instead of /3; for a constant barrier 
(B = 0), the first-crossing distribution is well-defined but the last- 
crossing distribution vanishes because there are continued cross- 
ings on all scales as 5 — > 5, . 

3 THE CORE MASS FUNCTION: RIGOROUS 
SOLUTIONS 

We have now derived the rigorous solution for the number of bound 
objects per interval in mass M, defined as the mass on the small- 
est scale on which they are self-gravitating. To apply this to a 
physical system, we need the collapse barrier B(S) and variance 
5 = a 2 (R) = a 2 (M). In HC08, B(5) is defined by the Jeans over- 
density pent (R) > ( c ? + v ?(B))/4-7rG, but the normalization of the 
background density, c„ and v, is essentially arbitrary. Moreover, 5 
is not derived, but a simple phenomenological model is used, and 
the authors avoid uncertainties related to this by dropping terms 
with a derivative in 5. In Hll, we show how S(R) and B(S) can 
be derived self-consistently on all scales for a galactic disk. For a 
given turbulent power spectrum, together with the assumption that 
the disk is marginally stable (Q = 1), S(R) can be calculated by 
integrating the contribution from the velocity variance on all scales 
R' >R: 



S(R) 



2 

"J : 



al{M[k])\W(k, R)\ 2 d\nk 



vj(k) 



'+K 2 k- 



(23) 
(24) 



where W is the window function for the density smoothing (for con- 
venience, we take this to be a Ar-space tophat inside k < l/R). This 
is motivated by and closely related to the correlation between Mach 
number and dispersion in turbulent box simulations (see |Padoan| 
et al.|f 



2008 



1997| |Passot & Vazquez-S emadeni 1998, Fede rrath et al. 
Price et al. 2011} . B(R) is properly given by 



B(R) = In 



Peril 

Po 



m 

2 



(25) 
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Figure 3. Variation in the predicted CMF/last-crossing MF with model as- 
sumptions. We compare the standard model from Fig.[T](p = 2, Aih = 30). 
Raising/lowering Aii, substantially slows/steepens the cutoff below the 
sonic mass. Assuming p = 5/3 (at fixed Mi,) also slows this cutoff. These 
variations have almost no effect on the high-mass slope, which is close to 
Salpeter in all cases. 



Pcrit 

Po 



Q_ 

2k 



1 + - 



h\ \o]{R) h 



.9 R 



Rj iaj(h) R 



where po is the mean midplane density of the disk, ft : 
for a constant-V c disk, and 

o 2 g (R)=c 2 + vl + {v 2 (R)) 



(26) 
-K/Q. = y/2 

(27) 



-1 



(28) 



The mapping between radius and mass is 

*w^4w [£+(i+£) «*(-?; 

It is easy to see that on small scales, these scalings reduce to the 
Jeans criterion for a combination of thermal (c s ), turbulent (v,), and 
magnetic (va) support, with M — (4n/3) pcrit/? 3 ; on large scales it 
becomes the Toomre criterion with M = 7rE cr i t i? 2 . 

Finally, we note that because the trajectories 8(x\R) defined 
above sample the Eulerian volume, the mass function is given by 

An PcritfM") „ ,„ x | d5 
dM 



M 



MM) 



dM 



(29) 



It is worth noting that, with S(R) and B(R) derived above, 
there are only two free parameters that together completely specify 
the model in dimensionless units. These are the spectral index p of 
the turbulent velocity spectrum, E(k) oc k~ p (usually in the narrow 
range p w 5/3 — 2), and its normalization, which we parameterize 
as the Mach number on large scales Ml = (v 2 (h))/(c 2 + v A ). Of 
course, we must specify the dimensional parameters h (or c s ) and 
po to give absolute units to the problem, but these simply rescale 
the predicted quantities. 

4 RESULTS 

Using these relations, we are now in a position to calculate the last- 
crossing mass function. Figure [T] shows the results of calculating 
the last-crossing distribution fi(S) for typical parameters p = 2 
(Burgers' turbulence, typical of highly super-sonic turbulence) and 
Mi, = 30 (typical for GMCs). First, we can confirm our analytic 
derivation in Eq.[l4] It is easy to calculate the last-crossing distri- 
bution by generating a Monte Carlo ensemble of trajectories S(R) 
in the standard manner of the excursion set formalism (beginning 
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at R — > oo and 5 = 0), and the details of this procedure are given 
in HI 1 ; here we simply record the last down-crossing for each tra- 
jectory and construct the mass function. The result is statistically 
identical to the exact solution. However, below the "turnover" in 
the last-crossing mass function, the Monte Carlo method becomes 
extremely expensive and quite noisy, because an extremely small 
fraction of the total galaxy volume is in low-mass protostellar cores 
(sampling <C 0.1M son i c requires ~ 10 10 trajectories). 

We compare the last-crossing mass function to the predicted 
first-crossing distribution - i.e. the mass function of bound objects 
defined on the largest scale on which they are self-gravitating. The 
two are strikingly different: they have different shapes, different 
power-law slopes, and the "characteristic masses" are separated by 
more than six orders of magnitude! Clearly, it is critical to rigor- 
ously address the "cloud in cloud" problem when attempting to de- 
fine either. 

Physically, in Hll we argue that the first-crossing distribu- 
tion should be associated with the mass function of GMCs, and 
show that it agrees very well with observations of the same. The 
last-crossing distribution, on the other hand, should correspond to 
the proto-stellar core MF - and by extension to the stellar IMF 
(although that conversion is uncertain and may well be mass- 
dependent, as discussed below). We see that here. As shown with 
the approximate derivation in HC08, the exact predicted mass func- 
tion reproduces all of the key observed features of the core/stellar 
MF. The high-mass behavior is an approximate power-law with a 
Salpeter-like slope d/V /AM oc M~ a with a slightly larger than 2; 
at lower masses the slope flattens and there is a subsequent low- 
mass turnover in good agreement with a Kroupa ( 2002 1 or Chabrier 
(2003) IMF. In Figure [2] we plot a direct comparison between the 
predicted CMF and observed starless/pre-stellar CMFs in differ- 
ent regions, which are less well-constrained than the stellar IMF 
but should correspond more closely with what we actually predict. 
The same behaviors are observed and agree well with our predic- 
tion without the need to invoke an uncertain rescaling factor (see 



e.g.|Stanke et al.|2006[ |Enoch et al.|2008l |Rathborne et al.|2009| 
Sadav oy"eTal.|2010||Be loche et al.|2011| >. 

Why do HC08 recover similar behavior, given that we have 
shown the first-crossing and last-crossing distributions can be so 
different? The generating equation for the mass function in HC08 
is the same as that of the original Press & Schechter ( 1974) deriva- 
tion of the halo MF; it neglects the cloud-in-cloud problem (treat- 
ing it as an ad hoc re-normalization). Modulo that renormaliza- 
tion, this is identical to the gi term in Eq. [14] This is a good ap- 
proximation when the integral/g2 term in Eq. [14] is small and/or 
proportional to fe(S). This occurs when dfi/dS 3> B(S)/S 3> 1; 
if the barrier decreases more rapidly than 5 (as R increases), 
then S(R) for a given trajectory will evolve slowly relative to 
B(R) and the probability of multiple crossings becomes small. In 
Eq. [14] the Poi term approaches a Dirac <5-function and the inte- 
gralterm^ -/,(S) + ©(dS/dS)" 1 , giving / £ (S) « (1/2) gi (5) + 
O (dB/dS) - 1 . At scales R/h<^l, relevant for the stellar IMF, S(R) 
is a weak function of R since most of the contributions to 5 come 
from the larger Mach numbers on larger scales, so this is satisfied^] 



2 This is very different from the reason that the original Pre ss & Schechter] 
(1974) mass function is correct to within a constant normalization, and from 
the nominal justification given in HC08 for assuming the same. That argu- 
ment relies on a constant barrier (dB/dS = 0) and applies only to the first- 
crossing distribution, for which the symmetry in P[q (broken in Pol) means 
that a trajectory at 8(S) = B(S) = Bq is equally likely to rise above or fall 
below the barrier. 



This approximation would break down severely at larger 
scales - in fact, the HC08 mass function, applied directly to the 
entire galaxy disk, would include a large part of the GMC mass 
function, and give the result that most of the mass in the core MF 
is in ~ 10 6 Mq "objects." But since the derivation in HC08 specifi- 
cally (and sensibly) truncates the problem before going to any large 
scales (focusing only on the regime R/h <C 1 and R/Rgmc < 1), 
this is avoided. Also, the normalization is different, but that is es- 
sentially treated as arbitrary by the authors. 

So over a limited regime, we can assume dB/dS 3> B(S)/S S> 
1, ignore the running of S(R), and approximate the mass function 
as 

dn _ pent (M) I din p cnt I (ln[p m t/po] +5 /2) 2 1 

dM ~ M 2 V2^> I dlnM I 6XP [ 25o 



(30) 



In the high-mass regime, the power-law slope is set by the 
run of p c ,it with R in the turbulence-dominated regime; for R in the 
range such that v< 3> c s and R <C h, S(M) ~ constant and p cnl oc 



R p 3 , M oc R p , so the approximate analytic scaling in Eq 
a slope 

3(l+p _1 ) , (3-p) 2 ln(M/M )-pln2 



30 



Q^highmass ' 



2S(M)p 2 



gives 



(31) 



where Mo ~ po h 3 ~ 10 s — 10 6 Mq for typical conditions. Based on 
the derivation in HC08, we would expect a somewhat "lognormal- 
like" turnover in the high-mass end of the CMF; however, this does 
not self-consistently calculate the running of S(M). In the exact cal- 
culation, the correction term oc 1 /S(M) is already small in the range 
of interest (« -0.1, for 5(M) ~ 10 and M/M ~ 10~ 6 - 10~ 4 ) 
and higher-order corrections nearly cancel, giving a very nearly 
Salpeter (a « 2.2 — 2.4) power-law behavior over ~ 2dex in mass. 

The turnover occurs below the sonic length, where thermal 
support both makes the equation of state "stiffer" and suppresses 
the running of the variance S(R), giving p cr j t oc R~ 2 , M oc R, so 



1.5+41n.M / ,-21n (M/M ) 



this 



S(M) 
becomes 



(32) 



CVlo 



For typical parameters 
0.161n(M/M somc ). 

Our solution also allows us to predict the normalization and 
mass scale of the mass function, without having to make any ad 
hoc assumptions to renormalize the problem about a given scale. 
The CMF begins to turn over below the sonic length, R s{ 



hM 



-2/0-1) 



with mass M sc 



(4V2~7r/3) .M 



sonic 



(2/3)cJfl so „ic/G w M (c/O.Skms" 1 ) 2 (fl sonic /0.1pc). The mass 
scale of the CMF the stellar IMF is therefore a natural consequence 
of Jeans collapse in a turbulent medium, with characteristic masses 
that depend only on the sound speed and sonic length (themselves 
related). This is a well-established result from both analytic work 
and numerical simulations (see e.g. Klessen & Burkert 2000, 2001 ; 
Bate & Bonnell 2005 1, and is true even if the gas is not isothermal 
(although our sub-sonic extrapolation is questionable in this case; 
see |Larsonj 1985 2005, Jappsen et al. 2005). There is, of course, 
some efficiency factor e cor e that relates the mass of a collapsing, 
bound protostellar "core" (which is what we actually predict) to 
the mass of a star. This may well depend on mass (and the thermal 
physics below the sonic length), which will introduce additional 
corrections between the IMF here and the stellar IMF (see §|5]be- 
low). However, the calculation here suggests that such corrections 
should be order unity, and arguments from observations and mod- 
els of outflows similarly suggest a modest e C ore ~ 0.5 jMatzner & 
|McKee|2000||Stanke et al.||20061|Alves et al.|[2007l|Enoch et al 
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|2008[ >. If constant (a significant assumption), this is comparable to 
the normalization uncertainty in Msonic from c s and i? son ic. 

In Fig. [3] we plot the predicted last-crossing distribution as 
a function of model parameters. We compare our "default" model 
(p — 2, Mi, = 30) to one with p = 5/3. For otherwise equal param- 
eters, this shifts M son i c to lower absolute values, but in units of the 
sonic mass the behavior is qualitatively the same. The high-mass 
slope is nearly identical - the difference predicted by the approxi- 
mate Eq.[30]or in HC08 is small to begin with, but is also largely 
canceled out by the second-order corrections from the proper last- 
crossing distribution. The low-mass turnover is much "slower," 
however, since the turbulent velocity declines more slowly relative 
to the sound speed (v, oc Z? 1 ' 3 instead of R 1 ^ 2 for p = 2), which 
enters into both the barrier and the run of S(R). Raising the Mach 
number on large scales, Mi,, has a similar effect, and slightly steep- 
ens the high-mass slope near the turnover mass. Decreasing the 
Mach number steepens the turnover below M SO nic, and for quite low 
Mi, also begins to manifest a high-mass cutoff. These are the dom- 
inant physical effects. Changing more subtle model choices - for 
example, using a Gaussian window function to smooth the density 
field instead of a i-space tophat, make relatively little difference. 

5 CAVEATS ON THE CMF-IMF CONVERSION 

What we have calculated here is the mass function of self- 
gravitating, non-fragmenting cores at a given "snapshot" or instant, 
for a galaxy with fully developed non-linear turbulence. It is im- 
portant to note that this is not yet the stellar IMF. As noted above, it 
is necessary to invoke some mean efficiency of ~ 30% in the con- 
version of core mass to stellar mass - presumably this comes from 
some combination of outflows and accretion histories, which could 
easily depend on mass; it is also possible that accreting protostars 
can be ejected from their environments, introducing a large core-to- 
stellar mass scatter. None of these physics are included in the cur- 
rent model; it may be possible to add them "on top" of the model 
here, but to do so it would be necessary to construct a more com- 
plex analytic model with some significant additional assumptions 
about how pre-stellar cores contract and grow, and how outflows 
affect the ambient medium. 

Although within the "snapshots" here, the cores have no sub- 
scales which are independently self-gravitating, there can certainly 
be successive fragmentation on smaller scales within the cores 
and/or stellar accretion disks as the cores contract and form stars 
(which will depend on the gas thermodynamics and stellar feed- 
back; see references above and e.g. Krum holz et al.|2009| [Peters] 
|et al.p OlO). Indeed this must occur, as a large fraction of stars 
are in binaries and high-mass stars are almost always members 
of higher order multiple stellar systems (Lada 2006; Zinnecker & 
Yorke 20071. If the multiplicity is mass-dependent (i.e. fragmenta- 
tion is not scale-free), this will also change the mass spectrum from 
CMF to IMF. There may be some prospects of including this in fu- 
ture, time-dependent excursion set models (discussed below), but 
if fragmentation is feedback-dependent or occurs within the stel- 
lar accretion disk, following it accurately depends on physics not 
included in the current model. 

Even if each core forms a fixed fraction of its mass into a 
fixed number of stars, the timescales for it to do so may be dif- 
ferent. What the IMF really samples is the relative "formation rate" 
of stars of a given mass; so if the timescale for collapse of massive 
cores is much longer than the timescale for collapse of low-mass 
cores, this can significantly change the shape of the IMF relative to 
the CMF (see the discussion in |Clark et al.|2007] >. At high masses 



M > M son i c , we do not expect the dependence of collapse timescale 
on mass to strongly change our results. For R > R son k and R<^ h, 

— 1/2 

the characteristic collapse time for cores scales as ~ r ayn oc p cr it oc 
R/vt oc R 1 / 2 oc M 1 / 4 . If the rate of star formation simply scales as 
oc Mcore/f collapse, this implies the bright end of the stellar IMF will be 
steeper than that of the CMF predicted here by M _1//4 . This is not 
a negligible correction, but it is also not very large (comparable to 
the difference between p — 5/3 and p = 2). Moreover, for p = 2 
it actually goes in the "right" direction: the CMF slope predicted 
(with second-order corrections, a ~ 2.15 — 2.2) is slightly more 
shallow than the canonical Salpeter slope, so this would bring the 
predicted IMF slope closer to observations. However, at low masses 
(<C M S onic), this effect is more important: ti yn oc R/c s oc R oc M, so 
the timescale correction can significantly modify the CMF-IMF re- 
lation. But in this regime, we already noted that our predictions are 
much more sensitive to the turbulent spectrum, and our simplifying 
assumption that the gas is isothermal, while introducing negligible 
corrections at R > i? SO nic, can have large effects on CMF in the sub- 
stellar (R < /? S onic) regime (see |Hennebelle & Chabrier|2009) . We 
therefore caution that while collapse times are critical to the CMF- 
IMF conversion at low masses, the gas thermodynamics is likely to 
be an equal or greater source of uncertainty (given the simplicity of 
our derivation here). 



6 DISCUSSION & CONCLUSIONS 

We have applied the excursion set model of the ISM, developed 
in Hll, to predict the mass spectrum of star-forming cores and, 
by (admittedly less certain) extrapolation, the stellar IMF. We note 
that the density field smoothed on a scale R about a random point 
x in the ISM, 5(x, R) may be self-gravitating on many scales. 
The largest of these scales defines the "first crossing distribution," 
which we show in HI 1 should be associated with the distribution of 
GMCs. If it is also self-gravitating on smaller scales, then this "par- 
ent" scale is not part of the core MF or stellar IMF, but will frag- 
ment into smaller units that collapse first: rather, the CMF should 
be associated with the "last crossing distribution," or the mass 
function defined on the smallest scales on which systems are self- 
gravitating. To our knowledge, this has not previously been studied 
in cosmological or ISM applications. This can be determined in a 
standard Monte Carlo excursion set approach, but we also derive 
the exact analytic solution for an arbitrary collapse threshold. We 
use this to predict the CMF, and find it agrees very well with that 
observed. The high-mass slope emerges as a consequence of tur- 
bulent support/fragmentation and is very close to a power-law over 
~ 2 dex, and very insensitive to the properties of turbulence for rea- 
sonable power spectra. The slope flattens at the sonic length/mass, 
~ c 2 R SO mc/G, as thermal pressure makes the effective equation of 
state more "stiff." The shape of this turnover is robust, but exactly 
how rapidly it occurs (in units of mass) depends on the Mach num- 
ber and turbulent spectral shape. 

Our calculation supports the conclusions of HC08, who de- 
rive the same behavior from a simpler "collapsed fraction" argu- 
ment, and do not attempt to calculate the variance in the density 
field or disk-scale effects. In many ways, the derivation here is re- 
lated to that work as that in Bo nd et aT7| ( |199l[ > is related to [Press] 
|& Schechter|(1974] >. We have provided a formal mathematical jus- 
tification, and exact solution, for the approximate mass function 
therein. Moreover we provide a means to calculate the normaliza- 
tion (of both space density and mass scale) and variance beginning 
at the parent disk scale, which removes many ambiguities. And 
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incorporating the argument into the full excursion set framework 
makes possible many additional calculations. 

We stress that the model used here to derive the CMF is exactly 
the same as that used in HI 1 to predict the properties of the GMC 
population and large-scale ISM structure. The only difference is 
that here we examine the last, instead of first, crossing distribution. 
Thus, beginning at the scale of an entire galactic disk, with the only 
"free parameter" being the turbulent spectrum, we derive first the 
distribution of GMCs, and simultaneously the MF of cores on much 
smaller scales near the sonic length, within those GMCs. This is, to 
our knowledge, the first analytic model that simultaneously predicts 
the CMF/IMF and large-scale ISM properties - both in remarkably 
good agreement with observations. 

There are a number of obvious extensions to our calculation. 
With a rigorous definition of the last-crossing distribution, it is 
straightforward to predict other properties derived for the GMC 
population in HI 1: the core size-mass and linewidth-size relations, 
the distribution of "voids" or "bubbles," the collapse rate of bound 
structures (and rate at which new mass is "supplied" into cores), 
and the dependence of these properties on turbulent spectral in- 
dex, Mach number, gas thermal properties, and more. With a well- 
defined fi (the last-crossing distribution; Eq. |14[ > is it possible to 
calculate the correlation function of individual stars and clustering 
properties of star formation. 

In HI 1, we also show how the excursion set formalism can be 
used to construct a Monte Carlo "merger/fragmentation tree," anal- 
ogous to extended Press-Schechter halo merger trees, with which 
it is possible to follow the time evolution of the ISM and GMC 
populations and build semi-analytic models for the same. Since the 
last-crossing is determined for all trajectories in the Monte Carlo 
formalism, it is straightforward to extend this to sample the IMF in 
a fully rigorous time-dependent manner, as a function of properties 
in different GMCs. With this formalism, GMCs, once bound, can be 
allowed to depart from the background flow, and this can in turn be 
allowed to change the CMF and ultimately the IMF. The IMF cal- 
culation could, in turn, be used to determine the per-GMC SFR and 
stellar mass distribution, which then informs whatever model of 
feedback is used to evolve the clouds themselves. One might imag- 
ine a model in which, for example, the GMC remains bound until a 
sufficient number of sub-regions have experienced last-crossing at 
a scale sufficiently large to form O-stars, which then dissociate the 
cloud. This also provides a means to generalize calculations such as 
that in Krumholz & McKee ( 2005 1 to rigorously derive the star for- 
mation rate, without having to assume an ad hoc collapse timescale 

It is also possible to follow time-dependent fragmentation. 
One concern with our "instantaneous" calculation (which essen- 
tially calculates last-crossing for the galaxy at a "snapshot" in time) 
is that a non-negligible amount of mass has its last crossing at scales 
much larger than any single star: see e.g. the overlap in the high- 
mass ft and low-mass tail of the first-crossing or GMC mass func- 
tion in Fig. [T] In a fully time-dependent model though, these re- 
gions would require at least ~ 1 dynamical time to contract, which 
amounts to many dynamical times for the denser sub-regions of the 
cloud. In this time, therefore, a large number of sub-regions might 
cross below the barrier - this is a purely analytic means to calculate 
time-dependent fragmentation of collapsing objects. This is likely 
to resolve the "problem" of these otherwise high-mass objects, but 
also to steepen the high-mass end of the predicted IMF relative to 
the CMF. Extension of our models to a fully time-dependent for- 
mulation may also provide a means to address the timescale cor- 
rections and fragmentation processes discussed in § [5] which are 



critical to understand the relation between CMF predicted here and 
stellar IMF. Including the effect of outflows is, however, non-trivial, 
and will require the introduction of fundamentally new physics in 
the models. 
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